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ABSTRACT 



We present a new measurement of the a-spectroscopic factor {So) and the 
asymptotic normahzation coefficient (ANC) for the 6.356 MeV 1/2"^ subthreshold 
state of ^''O through the ^•^C(^^B, ^Li)^^0 transfer reaction and we determine the 
a-width of this state. This is beheved to have a strong effect on the rate of the 
^^C(a,n)^^0 reaction, the main neutron source for slow neutron captures (the 
s-process) in asymptotic giant branch (AGB) stars. Based on the new width we 
derive the astrophysical S-factor and the stellar rate of the ^^C(a, nY^O reaction. 



At a temperature o: 



Caughlan fc Fowler 



■ 00 M K our rate is roughly two times larger than that by 



(119881 ) and two times smaller than that recommended by 
the NACRE compilation. We use the new rate and different rates available in 
the literature as input in simulations of AGB stars to study their influence on 
the abundances of selected s-process elements and isotopic ratios. There are no 
changes in the final results using the different rates for the ^'^C(q;, n)^^0 reaction 
when the ^^C burns completely in radiative conditions. When the ^'^C burns in 
convective conditions, as in stars of initial mass lower than ~2 Mq and in post- 
AGB stars, some changes are to be expected, e.g., of up to 25% for Pb in our 
models. These variations will have to be carefully analyzed when more accurate 
stellar mixing models and more precise observational constraints are available. 



Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: AGB and 
post-AGB 
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1. Introduction 

Approximately half of the elements heavier than iron in the universe are produced via 
a series of slow neutron capture reactions and competing /3-decays (the s-process). During 
the s-process, the neutron number density is relatively low, of the order of 10^ n cm~'^. 
When the flux reaches an unstable nucleus, it typically decays rather than capture another 
neutron and the 5-pro cess proceeds via the isotopes around the valley of /3-stability (e.g.. 



Burbidge et al. 



19571 ). The astrophysical sites of the g-p rocess are core He and shell C 



burning in massive stars for the elements lighter than Sr (jPignatari et al. 



20101), and the 



He intershell" of asy mptotic giant branch (AGB) stars for the elements between Sr and Bi 



( IGallino et al. 



19981 ). Stars with initial masses lower than roughly 9 Mq reach the AGB 
phase in the final phases of their evolution, when both H and He have been exhausted in the 
core leaving C and O in electron degenerate conditions. Production of nuclear energy occurs 
in the H and He shells, which are located between the core and the extended convective 
envelope and are separated by the thin He intershell layer. AGB stars experience thermal 
pulses (TPs) when the usually dormant He burning shell is suddenly activated. A large 
amount of energy is released, which drives convection in the He intershell. During TPs the 
star expands and cools and the H burning shell is inactive. While He burning turns from 
the convective to the radiative regime, and eventually switches off, the convective envelope 
can penetrate the underlying He intershell and carry to the surface the products of He 
burning, in particular carbon and the elements heavier than iron made by the s-process. 
This mixing process is known as the "third dredge-up" (TDU). After the TDU is ended, 
the star contracts and heats up again and H burning resumes until another TP occurs and 
the cycle is repeated. This sequence of events can occur from a few times to hundreds of 
times, depending on the stellar mass and the mass-loss rate. AGB stars suffer from very 
strong stellar winds, which erode the envelope roughly within a million years and shed the 
newly synthesized material mixed to the surface by the TDU into the interstellar medium. 
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Eventually, the C-0 degenerate core is left as a cooling white dwarf (see iHerwig II2005I . for 
a review on AGB stars). 



According to the current standard model ( Gallino et al. 



Goriely fc Mowlavi 



2000 



Lugaro et al. 



2003a : 



Crist alio et al. 



998; Busso et al. 1999 



2009al ). some protons must 



diffuse from the convective envelope into the He intershell at the end of each TDU in order 
to produce enough ^^C to account fo r the observed abun dances of the s-process elements 



at the surface of AGB stars (see also 



Busso et al. 



200ll ). A thin layer is then produced, 



known as the ^^C "pocket", which is rich in ^^C made via "'^^C(p, 7)^^N(/3+i/)-^'^C. When in 
this region the temperature reaches about 9 x 10^ K, the ^^C(a,n)^^0 reaction is activate d 
and generates neutrons tha t trigg er the s-process (jHoUowell fc Iben 



1988 



1988 



Gallino et al. 



Kappeler et al. 



1990 



201l|). 



section ( 


Sekharan et al. 


1967; 


Davids et al. 


1968; 


Bair et a 




1973; 


Kelloee et al. 


Drotleff et al. 


1993; 


Brune et al. 


1993; 


Harissopulos et al. 


2005; 


Heil et al. 


2008) 



1989 



These measurements have been performed at energies down to 270 keV, whereas the Gamow 
window is at 190 ± 40 keV, corresponding to a temperature of 100 MK. Since this energy 
is far below the Coulomb barrier, the reaction cross section is extremely small and direct 
measurement is sensitively limited by background signals and very difficult to perform 
in laboratories on the Earth's surf ace. While a measurem ent has been proposed at the 



underground laboratory of LUNA ( ICostantini et al. 



20091 ). at present, the experimental 



cross sections have to be extrapolated below 2 70 keV. A microscopic cluster model analysis 



of the ^^C(Q;,n)^^0 and ^^C{a,a) reactions by iDescouvemont I (119871 ) suggested that this 
extrapolation is critically affected by the 1/2"*" subthreshold resonance in ^^O {Er^ = 6.356 
MeV, just 3 keV below the a-decay threshold). The contribution from this resonance 
depends strongly on the a-width of the 1/2"*" state in ^^O, which can be derived from the 



- 5 - 



spectroscopic factor [So) or the asymptotic normalization coefficient (ANC) of a-cluster in 
this state. 

The Sa and the ANC can be determined from the angular distribution of the direct 
a-transfer reaction using distorted wave Born approximation (DWBA) or coupled reaction 
channels (CRC) analysis. Although th ree indirect measurern e nts via the (^Li , d) or the 



(^Li, t) system have be en performed by 



Pellegriti et al. 



Kubono et al. 



fl2003h . 



Johnson et al. 



(120061), and 



(|2008[ ) to study the Sa or the ANC of the 1/2+ state, there still exists a 
significant discrepancy of up to a factor of ~30 in the derived Sa and ANC. Therefore, it 
is interesting to perform a new measurement of the S^ and the ANC via an independent 
transfer reaction. In addition, it is necessary to understand the impact of the different 
resulting ^^C(a;,n)^^0 rates on the s-process nucleosynthesis in AGE stars. 

In this paper we determine a new stellar rate of the ^^C(q;, ?t,)^^0 reaction and 
incorporate it in calculations of the s-process nucleosynthesis in AGE stars. First, we 
measure the angular distribution of the ^^C(^^E, ''Li)^'^0 reaction to determine the Sa and 
the ANC for the 1/2+ state in ^"^O. Using this experimental ANC we derive the a-width 
for the 1/2+ subthreshold resonance, which is currently the most uncertain parameter for 
determining the ^'^C(a, n)^^0 rate. Finally, we use the new rate and different rates available 
in the literature as input for simulations of AGE stars to study their influence on the 
abundance of some selected s-process elements and isotopic ratios. 



2. Measurement and evaluation of the ANC 

2.1. Experimental procedure 

The measurement of the angular distribution for the ^^C(^^E, ''Li)^^0 reaction was 
performed at the HI- 13 tandem accelerator of China Institute of Atomic Energy (CIAE) 
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in Beijing. We used a ^^B beam with an energy of 50 MeV leading to the production of 
the excited states in ^^O at = 3.055 MeV, 3.843 MeV, 4.554 MeV and 6.356 MeV. The 
angular distribution of the elastic scattering for the entrance channel (^^B + ^^C) was also 
measured. We used a self-supporting ^^C target with a thickness of 75 ±6 jig/cm^ and 
an initial purity of 88%. A 26 MeV ^Li beam was also dehvered for measurement of the 
exit channel (^Li + ^^0) elastic scattering. Natural silicon monoxide of 86 ±7 fig/cvo? was 
prepared onto a 40 ±3 //g/cm^ carbon foil, serving as the oxygen target. In addition, a 
self-supporting ^^C target of 66 ±5 jig/cw? was used for cahbration of the focal plane and 
background subtraction during the whole experiment. 

To monitor the possible buildup of ^^C, the ^^B elastic scattering on the ^^C target was 
measured at the start and at the end of the measurement for each angle. This showed that 
the ^^C buildup was negligible compared to the initial impurity in the ^^C target, possibly 
due to the rather low gas pressure (< 10~^ mb) in the reaction chamber. To determine 
the absolute amount of ^^C in the ^^C target we measured the angular distributions of the 
iig _|_ 12q elastic scattering with both the natural ^^C target and the enriched ^^C target. 
The absolute amounts of ^^C and ^^C in the ^^C target were found to be 9.0 ±0.7 and 
66 ±5 //g/cm^, respectively. 

The beam current was measured by a Faraday cup covering an angular range of ± 6° 
in a laboratory frame and used for the absolute normalization of the cross sections at 
^lab > 6°. The Faraday cup was removed when measuring the cross sections at ^lab < 6°. 
A Si AE" — E telescope located at ^lab — 25° was employed for the relative normalization 
of the cross sections at ^lab < 6° by measuring the elastic scattering of the incident ions on 
the targets. In addition, the ratio of current integration in the Faraday cup to the elastic 
scattering events was measured at the start and the end of the measurement for each angle 
^lab ^ 6° by restoring the Faraday cup. The ratios changed by less than 2%, which thus 
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led to a reliable normalization of the cross sections at Oi^h < 6°. 

The reaction products were focused and separated by a Q3D magnetic spectrograph 
and recorded by a two-dimensional position-sensitive silicon detector (PSSD, 50 x 50 mm) 
fixed at the focal plane of the spectrograph. The two-dimensional position information 
from the PSSD enabled the products emitted into the acceptable solid angle to be recorded 
completely. The energy information from the PSSD was used to remove the impurities with 
the same magnetic rigidity. 

Due to the presence of ^^C in the ^^C target, the ^Li events from the ^^C(^^B, ^Li)^^0* 
(6.356 MeV) reaction were mixed with those from the ^'^C(^^B,'^UY^O* (6.917 MeV) 
reaction. To evaluate this background the (^^B,^Li) reactions were measured for both the 
^^C and ^^C targets at each angle with the same experimental setup. As an example, in 
Fig. [T]we display the focal-plane position spectra of ^Li at ^lab = 10° from the (^^B,^Li) 
transfer reactions. The background from ^^C accounts for approximately half of total events 
from the ^'^C target. After background subtraction and beam normalization, the angular 
distributions of the elastic scattering and the ^'^C(^^B, ''Li)^'^0* (6.356 MeV) reaction were 
obtained, as presented in Figs. |2] and |3l In Fig. [3] we also display the angular distributions 
of the i3C("B,'^Li)i'^0 reaction leading to the other three states {E^ = 3.055 MeV, 3.843 
MeV, and 4.554 MeV). Note that the measurements for these three states were not affected 
by the ^^C background. 



2.2. Extraction of the ANC 



The finite-r ange distorted wa ve Born approximation (DWBA) method with the 
FRESCO code iThompson 1988 ) was used to analyze the experimental angular 
distributions of the transfer reaction. The optical model potential (OMP) parameters 
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for the entrance and exit channels were obtained by fitting the experimental angular 
distributions of the ^^B + ^^C and ''Li + ^^0 elastic scattering, respectively (Fig. Full 
complex remnant term interactions were included in the transfer reaction calculations. The 
parameters of the core-core (''Li + ^^C) potential of 



Cook et al. 



(119871 ) gave reasonable 



account of the elastic scattering of "^Li from ^^C at 34 MeV. The parameters used in the 
DWBA calculations are listed in Table [TJ 

To obtain the spectroscopic factor {Sa) and the ANC of a-cluster in ^^O the 
spectroscopic amplitudes of a-cluster in the ground state of ^^B need to be fixed. The 
single-particle wave function describing the relative motion between a-cluster and ''Li core 
in the ^^B ground state can have two components denoted by quantum numbers NLj = SSq 
and 2D2, respectively, where N is the number of radial nodes that include the origin but 
not the infinity and L and j are the orbital and total angular momenta, respectively. The 
spectroscopic amplitudes of these two co mponents are —0.509 and 0.629, respectively, from 
a shell model calculation (IKurath 1 119731) . and —0.638 and —0.422, respectively, from the 



translationally invariant shell model fiRudchik et al. 



20051 ). In the present analysis, both 



sets of spectroscopic amplitudes were used and the resulting difference was incorporated in 
the total uncertainty of our result. 

The a-cluster single-particle wave functions were calculated using conventional 
Woods-Saxon potentials whose depths were adjusted to reproduce the binding energies of 
a-cluster in the ground state of ^^B and in the four states of ^''O. The quantum numbers 
NLj for a-cluster in the 6.356 MeV 1/2+ state of ^^O were fixed to be 4Pi by the oscillatory 
energy conservation relation 2(A^ — 1) + L = Ylt=i'^i''^i ~ 1) + h, where nj,/j are the 
corresponding single-nucleon shell quantum numbers. In addition, the quantum numbers 
for a-cluster in the other three states, 3.055 MeV 1/2", 3.843 MeV 5/2~, and 4.554 MeV 
3/2~, were determined to be 45*0, SD2, and 3D2, respectively. 
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The geometry parameters, radius ro and diffuseness a, of the Woods-Saxon potential for 
a-cluster in ^^B were adjusted to give the root-mean-square (RMS) radius (^P) = 3.204 
fm) of the a-cluster wave function. This was calculated using the following relation between 
cluster sizes and their mean distances within a nucleus, 



m 



B 



m 



B 



"^He"^Li^^2^^ 



m 



B 



where the RMS radii of ^He, ''Li, and ^^B were taken to be 1.47 fm, 2.384 fm, and 2.605 fm. 



respectively (ILiatard et al. 



19901). The resulting parameters are tq = 0.92 fm and a = 0.65 
fm. We investigated the dependence of the calculated Sa and the ANC on the geometry 
parameters for ^^B. With a diffuseness between 0.65 — 0.75 fm, the radius was adjusted to 
reproduce the RMS radius of 3.204 fm. The impact of this change on the Sa was found to be 
less than 2%. We also investigated the dependence of the and the ANC on the geometry 
parameters for ^^O using the same range of ro (0.9 — 1.06 fm) a nd a (0.60 — 0.76 f m) th at 



Pellegriti et al. 



( 120081 ). 



was selected by maximum likelihood function set at 3a level by 
Using steps of 0.02 fm for both tq and a, 81 sets of geometry parameters were obtained and 
used to calculate 81 values of the Sa and the ANC. Their standard deviations were taken as 
the uncertainty deriving from the geometry parameters for ^^O. The parameters ro = 1.00 
fm and a = 0.76 fm provide the best description for the angular distributions of the four 
states (see also Fig. [3]). 

In Fig. [3] we display the calculated angular distributions normalized to the experimental 
data for the ^'^C(^^B,'^Uy'^0 reaction populating the 3.055 MeV, 3.843 MeV, 4.554 MeV, 
and 6.356 MeV states in ^''O. For the 3.055 MeV state the Sa factor was found to be 
0.1 9 ±0.06, compatible with the values obtained from the (^Li, d) reaction (Sg ^ 0.18 — 0.3) 
hoO^j ] and from the (^Li,t) reaction (^^ = 0.27 ± 0.05) by 



Keelev et al. 



Pellegriti et al. 



(120081 ). For the 3.843 MeV and 4.554 MeV states the 5*^, factors were found to be 



0.078 ±0.025 and 0.060 ±0.019, respectively. These disagree with the values of 0.19-0.34 
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and 0.27 — 0.48, respectively for the two states, given by lKeeley et al. I (|2003[ ). who adopted 
geometry parameters tq = 1.25 fm and a = 0.65 fm for ^^O and used a coupled-channel 
calculation to fit the data. When using these geometry parameters we were not able 
to reproduce our measured ^•^C(^^B, ^Li)^^0 angular distributions of the four states and 
obtained an extremely small value of ma ximum likelihood function. For the 4.554 MeV 



Pellegriti et al. 



(120081 ) ■ who found an Sa factor 



state our result is consistent with that of 
of 0.10 ±0.05 via the (^Li,t) reaction. 

The Sa for the 6.356 MeV 1/2+ state of ^^O was derived to be 0.37 ±0.12. The 
error results from the statistics (23%), the target thickness (8%), the uncertainties 
from the spectroscopic amplitudes (3%) for ^^B, the geometry parameters (2%) for 
^^B, and the geometry parameters (20%) for the 1/2+ state of ^^O. The square of the 
Coulomb modified ANC (C^) was then extracted to be 4.0 ±1.1 fm~^ using the relation, 
C*^ = SaR'^ (p{R)'^ /W {RY , where (p{R) is the radial single-particle wave function for 
a-cluster in the 1/2+ state of ^'^O, and W{R) = W{R)T{L ± 1 ± r/) is the Coulomb modified 
Whittaker function, with r(L± 1 ±77) being gamma function and rj the Coulomb parameter. 
The uncertainty in the ANC of 27.5% is smaller than that in the Sa of 32.4%. This is 
because the variation of the geometry parameters yields a change in (j){R) that is opposite 
in sign to the change in Sa so the uncertainty from the ^^O geometry parameters in of 
12% is smaller than that in Sa of 20%. 



2.3. Different evaluations of the ANC 



Three independent measurements in addition to the present work have been performed 
to date to study the ANC or the S^ of the 1/2+ state in -'^'^O. A very small spectroscopic 
factor (5*0- ~ 0.011) was found via me asurement of the ^^Cf ^Li, (i)^^O angular distribution 



with an incident energy of 60 MeV (iKubono et al. 



2OO3I ). This indicated that the 



- 11 - 



contribution of the 1/2+ subthreshold resonance is neghgible . However, a re a nalys is of the 



same expe rimental data showe d that the DWBA analysis of 



Kubono et al. 



be flawed (IKeeley et al. 



fl2003[ ) could 



20031 ). These authors derived larger Sa factors of 0.36 and 0.40 



Johnson et al. 



via DW BA and coupled reaction channels (CRC) calculations, respectively. 
( 120061 ) measured the ^Li(^^C, (i)^''0 angular dist ribution at sub-Coulona b energies (8 and 



Pellegriti et al. 



(120081) measured the 



8.5 MeV) and derived = 0.89 ±0.23 fm-^ 
^^C(^Li, t)^^0 angular distribution with two incident energies (28 and 34 MeV). The Sa and 
were found to be 0.29 ±0.11 and 4.5 ±2.2 fm~^, respectively. 

In Table [2] we list the and C*^ values from the present study and from the 
literature sources me n tioned above. G o od ag reement for the Sa and is found between 



Keeley et al. 



mm 



Pellegriti et al. 



(120081 ). and the present work, which used three 



different transfer systems and covered an energy range of 4 — 10 MeV/u. This indicates 
that the Sa and ANC for the 1/2"^ state of ^^O should not significantly depend on the 
selection of transfer systems and incident energies. On t he other hand, the obtained 
from the ^Li(^^C, d) reaction at sub-Coulomb energies by 



Johnson et al. 



(I2006h is about 



Pellegriti et al. 



(120081). 



five times smaller than that obtained in the present work and by 
The sub-Coulomb (^Li, d) a-transfer cross section was recently remeasured to understand 
the source of this discrepancy and the data analysis is in progress (G. V. Rogachev 2012, 
private communication). It remains to be seen if this analysis will result in a revised value 
of 
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3. The ^^C{a,ny^O reaction rate 



The astrophysical S-factor of the ^^C{a,nY^O reaction via the resonances can be 
calculated with the Breit-Wigner formula, 



S{E) = 71 



2Jr + 1 



r,(E)r„(E + Q) 



■ exp(27r?7), 



(2) 



2/i (2Jp + l)(2Ji + 1) - Enf + iVtot/2f 

where /i is the reduced mass of the a + ^'^C system; Er represents the resonance energy; Jr. 
Jp, and Jt are the spins of the excited states in ^^O, a, and ^^C, respectively; Fq,, r„, and 
Ttot denote the a-, neutron-, and total widths, respectively; and Q is the reaction Q- value 
of ^^C{a,ny^O. 

For the resonances above the a-threshold, the energy dependence of these partial 
widths is given by 

PkiE) 



T^{E) = T^{E 



R) 



and 



r„(E + g) = r„(E^)( 



E + Q \ '/+1/2 



(3) 



^Er + Q 

where Ta{ER) and r„(i?/j) denote the experimental a- and neutron-widths, respectively 



(4) 



PiXE) is the a-penetrability; and li and If are the orbita 
neutron in the excited states of ^^O, respectively (see, e.g. 



angular mo menta for a and 



lliadis 



20071 ). 



For the subthreshold resonance, the energy dependence of the neutron-width can be 
also obtained by Eq. Hj while the dependence of the a-width is expressed as 



T^{E) = 2^lPiXE). 



(5) 



Here the reduced a-width 7^ can be given by 



7o 



Sa4>{Rc 



2/xi?, 



2 

cj 5 



(6) 
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and was extracted to be 12.7 ±3.5 keV at the channel radius Rc = 7.5 fm. This large 
radius was chosen to reach the Coulomb asympto tic behavior of 4>{R), as suggested by 



Descouvemont 



fll987f ) and 



Pellegriti et al. 



fl2008h . 



In Table [3] we list the resonant parameters for the 1/2"*" subthreshold state employed 
to obtain the astrophysical S-factor of the 1/2"^ subthreshold resonance in ^^C(a,n)^^0. 
The uncertainties in the S-factor were investigated by varying 7^, Er, and r„(i?ij). We 
also investigated the dependence of the S-factor on the channel radius by changing Rc 
from 6 fm to 7.5 fm. The S-factor at the Gamow peak of 190 keV, ^(IQO), was derived 
to be (8.4 ± 2.3) X 10^ MeV b. The error is dominated by the uncertainty in the resonant 
parameters for the 1/2"*" subthreshold state. The uncertainty from the channel radius is 4% 
in S'(190), varying slightly with the energy. 



We calculated the total S-factor of C{a,n) O 



17 



O states up to 8.342 MeV from the compilation of 



by including the p roperties of the 



Tillev et al. 



(1l993l ). and considering 



their interferences. However, the resulting tota. 



from the direct measurements of 



Drotleff et al. 



the lowest energies to date, and of 



HeiletaL 



S-fact or did not agree with the data 



19931) . who measured the S-factor at 



(120081 ). who performed a systematical 



experimental verification of neutron efficiency over a range of well-defined energies. 
Moreover, the resonant structure near the 3/2"*" resonan c e at ■E b = 0.842 MeV was not 



well reproduced using the parameters from 



Tilley et al. 



(I1993I ). Hence, we adjusted the 



a- and neutron-widths for this state to provide the best fitting of the experimental data. 
The fitting resulted in broa der partial widths i TniEfj) = 0.08 keV, Tn{Eji) = 340 keV] than 



the recommend ed values of 
agreement with iPellegriti et al. 



Tillev et a. 



(I2OO8I) and 



fll993h \TJEr,) = 0.07 keV, T„iEn) = 280 keV], in 



HeiletaL 



(120081 ). All resonant parameters 



used in the present calculations are listed in Table [3l 



In Fig. Hlwe display the resulting astrophysical ^^C(a,n)^^0 S-factor as a function 
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of the energy in the center of mass frame [Ec.m.)- The uncertainty mainly results from 
the resonant parameters of the 1/2+ subthreshold state since the resonant parameters of 
the states above the a-emission threshold are relatively well constrained by fitting the 
experimental data. The contribution of the 1/2"*" resonance results in a clear increase of the 
S-factor at lower energies. At the Gamow peak of 190 keV the 1/2"'" subthreshold resonance 
dominates the ^^C(a, n)^^0 reaction. At this energy the contribution of this resonance 
accounts for 61% of total S-factor. 

The astrophysical ^^C{a,ny^O reaction rate was calculated with the present S-factor 
using: 



NA{av)=NA(— 



\ 1/2 



S{E) exp 



E 



dE, 



(7) 



Vvr/i/ (A;T)3/2 

where is Avogadro's number. In Table H] we list the present adopted rate as function of 
the temperature together with its upper and lo wer limits. In Fig. [5] we co mpare the present 



rate with the previo us compilations by CF88 ( jCaughlan &: Fowler 
1 1 I 



1988h and NAC RE 



(|Angulo et al. 



Kubono et al. 



19991) and with other rates available in the 



2003 



Johnson et al. 



2006 



Pellegriti et al. 



iterature (jPenker et a. 



2008 



Hei 



Denker et al. 



et al. 



1995 



20081) at the 



(119951 ) and it is about 



temperature of 100 MK. Our rate agrees with that of 

two times smaller than the NACRE recommended rate, but within its lower limit. Both 
the rates by NACRE and 



Denker et al. 



(119951 ) were derived from an extra polation of the 



lowest energy experimental data available to date of 



Drotleff et al. 



( 119981 ). There is a 



large discrepaii c y of u p to a factor of 2 between t he present rate and the rates by CF88, 



Kubono et al. 



(l2003h . and 



Johnson et al. 



toom due to the fact th at CF88 did not take 



into account the con t ributi on of the 1/2"*" subthreshold resonance and 



Kubono et al. 



and 



Johnson et al. 



(120031) 



(120061 ) found a significantly smaller contribut ion than ours o 



" this 



resonance to the rate. Our rate is in good agreement with that of 



Pellegriti et al. 



mm, 



who found the contribution of the subthreshold resonance similar to ours, and is consistent 
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with that of iHeil et aL I ( 120081 ). who used an extensive multichannel i?- matrix analysis to 
constrain the possible contribution from subthreshold resonances by taking into account all 
open reaction channels for the ^^C + a and + n systems. 

For conveniency, we fitted our rate with the expression used in the astrophysical 



reaction rate library REACLIB ( iThielemann et al. 



1987 



Rauscher fc Thielemann 



20011 ) 



Na{ctv) = exp[ai + a2Tg ^ + asTg + a^Tg^^ + a^Tg + a^Tg^^ + a7ln{Tg)] 

+ exp[as + agTg"^ + aioTg^^^ + auTg^^ + a^Tg + ai^Tg^^ + ln(T9)]. (8) 

Here Tg is the temperature in units of 1 GK. The coefficients for our adopted value and 
lower and upper limits of the ^^C{a,ny^O reaction rate are listed in Table HI The overall 
fitting errors are less than 7% at temperatures from 0.04 to 10 GK. 



4. Astrophysical implications for s-process nucleosynthesis in AGB stars 



To assess and understand the impact of the new rates of the ^^C(a, ?7,)^^0 reaction 
we used a post-processing code where a nuclear network of 320 species (from H to Bi) 
and 2,336 reactions is solved simultaneously with mixing within the star when convective 
regions are present. We used stellar structure inputs, such as temperature, de nsity, and 



convective velocity, calculated previously by the Stromlo stellar struc ture code (jLattanzio 



Vassiliadis fc Wood 



9861) including mass-loss on the AGB phase with the prescription of 
(119931 ). We included the formation of the ^^C pocket by artificially allowing an exponentially 
decreasing proton profile to form just below the base of the convective envelope at the 
end of each TDU episode over a mass of 0.002 Mq (roughly 1/10**^ of the mass of the He 



intershell, in the low- mass AGB models considered here). The details of this 
the codes used to compute the models have been described previously (e.g.. 



2004 



Karakas et al. 



20091) 



procedure and 



Lugaro et al. 
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We considered four models, which are summ arized in Table [5l a M = 3 Mq model of 



metallicity Z =0.02 si milar to that discussed 



model of Z =0.01 from 
oiZ = 0.0001 from 



Karakas et al. 



Lugaro et al. 



by 



Lugaro et al. 



fl2003af ). a M =1.8 Mq 



fl2010f ). and M =1 Mq and M = 1.5 Mq models 



(120121 ). We selected these models for three reasons. (1) 
Models of higher masses probably do not experience the ^^C neutron source as protons 
burn while being inges t ed in the He i ntershell, preventing the formation of the ^^C pocket 



( iGoriely fc Siess 



2004 



Herwig 



20041 ). It is still not clear what the effect of this burning is 



on the whole stellar structure and the s-process, but it seems reasonable to us for the time 
being to not include a ^^C pocket for masses larger than ~4 Mq. In these massive AGB 
stars, the temperature in the TPs exceeds 300 MK so that the ^^Nefa, n^^^Mg reacti on is 



activated and likely plays the role of the main neutron source (Ivan Raai et al. 



The C in the pocket normally burns before the onset of the following TP ( 



2012). (2) 



Straniero et al. 



19951 ). however, in some cases the temperature may not be high enough for this to happ en 



and so me C could be left over to burn in the following TP (see also 



Cristallo et al. 



2009al ). In our lowest-mass models the first few ^'^C pockets are ingested in the following TP 
while they still contain a relatively large amount of ^^C. These models are thus qualitatively 
different from the M = 3 Mq model where ^^C always burns before the onset of the next 
TP. In details, we have 6, 15, and 2 TDU episodes in the M = 1.8 Mq model of Z =0.01, 
and the M = 1.5 Mq and M = 1 Mq models of Z =0.0001, respectively. These numbers 
represent also the number of ^^C pockets introduced in each models, since a ^^C pocket is 
introduced after each TDU episode. Of these ^^C pockets, the first 4 and 3 for the M = 1.8 
Mq and M =1.5 Mq models, respectively, are ingested in the following TP while the 
abundance of ^'^C by number is still higher than 10~^. The same occurs for both the ^^C 
pockets introduced in the M = 1 Mq model. As the AGB star evolves the temperature in 
the He intershell increases so that this effect disappears for later ^'^C pockets. The exact 
number of ingested ^'^C pockets relative to their total number increases with decreasing the 
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stellar mass and increasing the metallicity. This is because the temperature decreases for 
lower masses and higher metaUicities and thus the chance of some unburnt ^•^C surviving to 
be ingested in the next TP is higher. (3) Finally the low-mass and low-metallicity models (1 
Mq and 1.5 Mq, Z = 0.0001) present a further regime for the activation of the ^^C neutron 
source. Proton ingestion episodes occur in the first few TPs during which a small amount of 
protons are ingested directly ins ide the TPs and thus some extra ^'^C is prod uced and burnt 



inside these convective regions (I Crist alio et al 



2009b 



Lugaro et al. 



2OI2I). In Fig. Iwe 



schematically illustrate the location in time and space where the different regimes operate 
within the thermally-pulsating structure of an AGB star and indicate the models where the 
regimes occur. 

In Fig. [7] we display the variations of the abundance of some selected s-process 
elements and isotopic ratios at stellar surface at the end of the evolution obtained from 
our models using the present ^^C(a;,n)^®0 rate with respect to the two compilations, 



CF88 and NACRE, respectively (see 



Denker et al. 


( 


1995] 


from 


Johnson et al. 



ig. [5] for a cornpariso n of reaction rates 



Pellegriti et al. 



mm 



Heil et al. 



The 



mm . and 



1995) agree with the present rate within 10% for all the models. The rate 



(I2OO6I ) produces larger variations , but always wi t 



lin t 



lose reported 



for the CF88 and NACRE rates. Using the rate from iKubono et al. I (120031 ) on the other 
hand produces variations outside of the ranges in Fig. [TJ e.g., Pb is up to 40% lower than 
the value obtained using the present rate. We chose to plot the elements representing the 
first s-process peak (Sr), the second s-process peak (Ba), and the third s-process peak 
(Pb). These can be observed in stars and provide a description of the overall s-process 
distribution, which is a function of the total time-integrated neutron flux. We further plot 
two isotopic ratios: ^®Zr/^^Zr and ^^Kr/^^Kr. These are sensitive to the activation of the 
branching points at ^^Zr and *^Kr, respectively, and thus to the local value of the neutron 
density. These ratios can be measured in meteoritic Stardust silicon carbide (SiC) grains 
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that originated in C-rich AGB stars. We also report on the production of fluorine, however, 
we find that its production is not significantly sensiti ye (by ^10% at mos t) to the choice of 



the C{a,n) O reaction in all four models (see also 



Lugaro et al. 



20041 ) 



In the 1.8 Mq, Z =0.01 model we find that the overall s-process distribution changes 
when different ^^C(a,n)^^0 rates are used. The present rate produces 16% more Ba and 
25% more Pb than the slower CF88 rate. On the other hand, the present rate results in a 
14% lower abundance of Pb than the faster NACRE rate. These differences are a signature 
of a different efficiency of the s-process. A lower efficiency produces less Ba and Pb than 
Sr because for the neutron capture fiux to reach Ba and Pb higher number of neutrons 
per Fe seed are required. Slower rate results in more ^^C left over in the pocket to be 
engulfed in the following TP and the s-process is less efficient when ^^C is engulfed in the 
TP than when ^^C burns radiatively for two reasons, (i) In the radiative ^'^C pocket there 
is no mixing among the different layers of the pocket and thus ^^N nuclei produced in the 
region of the pocket where the initial number of proton is higher than 0.01 are not mixed to 
the ^^C-rich layer. This produces the highest s-process efficiency because ^^N is a neutron 
poison via the ^^N(n, p)^^C reaction, which removes neutrons from being captured by Fe 
seeds and their progeny. In the convective TP instead ^^C nuclei are mixed with the ^^N 
present in the pocket as well as in the H burning ashes, (ii) The neutrons in the radiative 
^^C pocket are released over a very small mass (0.002 Mq in our models), there are more 
neutrons present locally per initial Fe seed and the fiux of neutron captures can reach to 
the heaviest elements up to Pb. In the convective TP instead the neutrons are released over 
a larger mass (~ 0.01—0.02 Mq in our models), there are less neutrons per Fe seed resulting 
in the production of the lighter s-process elements like Sr. The isotopic ratios plotted in 
Fig. [7] behave in the opposite way of Pb because the local neutron density is higher, and 
thus branching points are more activated, when ^^C is ingested in the TPs due to the higher 
temperature resulting in a shorter burning timescale. 



- 19 - 



Interestingly, the low-metallicity Z — 0.0001 models present the opposite result to the 
Z = 0.01 model. E.g., in the 1 Mq model the present rate produces 25% less Pb than the 
slower CF88 rate. These stellar models have lower masses and a much larger fraction of 
^^C burns while ingested in the TP, together with ^^N, due to both incomplete burning of 
the ^^C pocket before the next TP and the proton ingestion episodes. As mentioned above, 
^''N acts a neutron poison via the (n, p) channel, but at the same time it is destroyed 
by a-captures in the TP. Hence, the timescale of ^^C ingestion and burning with respect 
to that of ^^N has an impact on the final neutron flux. This could explain why slower 
^^C(Q;,n)^^0 rates produce more free neutrons since if neutrons are released at later times 
there is less ^^N to capture them. In these conditions, one should also consider the effect of 
recycling of the protons made by the ^^N(n, p)^^C reaction. These protons can be captured 
by the abundant ^^C to make more ^^C, but also by ^^C itself. It is difficult to evaluate 
analytically the final outcome of all these combined effects. Overall, our numerical models, 
which solve simultaneously mixing and burning in the TPs, indicate that the present rate 
results in a lower overall neutron flux than the CF88 rate. This is accompanied by a higher 
or a lower neutron density, depending on the stellar mass and the exact ratio of ^^C burning 
convectively to ^^C burning radiatively. 17% lower and 22% higher ^^Kr/^^Kr ratios, which 
are extremely sensitive to even small changes in the neutron density, are obtained for the 1 
Mq and 1.5 Mq models, respectively. 

In the 3 M©, Z — 0.02 model we do not find any variations in the resulting abundances 
within 5%. This is because in this model ^^C always completely burns before the onset 
of the next TP. The total integrated neutron flux is thus determined only by the initial 
amount of ^^C in the pocket and not by how fast it burns. We expect the same behavior for 
more massive AGB model if a ^^C pocket was to be considered in these cases. The isotopic 
ratios in this model are also not sensitive to the ^^C(a,n)^^0 rate as they depend mainly 
on the neutrons released in the TPs when the ^^Ne neutron source is marginally activated. 
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It is interesting to compare the predicted abundances of the s-process elements and 
isotopic ratios with the observations. We find that some differences up to 25% are predicted 
by the s-process models, particularly for Pb, when varying the ^^C burning rates. Pb can be 
measured in carbon-enhanced metal-poor (CEMP) stars in the Milky Way halo, which are 
believed to carry the chemical signature of mass transfer from a more massive companion 
during its AGB phase. However, the spectroscopic abundances in s-process-enhanced 



stars are current 



Masseron et al. 



y de termined with typical uncertainties around 80% (see Table 2 of 



20101), larger than the differences found here. In addition, the formation 
of the ^^C pocket itself is very uncertain and there is no agreement on exactly which 
mechanism drives it, and on the impact of processes such as rotation and magnetic fields. 
Also the details of the proton ingestion episodes found in our low-metallicity models depend 
on the treatment of convective borders in stars, one of the largest uncertainties in stellar 
modeling. Future observations with high resolution and high signal to noise ratio spectra of 
low-metallicity stars, and further understanding of the formation of the ^^C pocket and the 
proton ingestion episodes will be all important for better constraining the s-process model 
predictions. 

Isotopic ratios affected by branching points are indicative of the ^^C burning rate, even 
though also in this case, other uncertainties including the neutron-capture cross section 
of the unstable nucleus at the branching point may have a more important effect. The 
^^Zr/^'^Zr ratio has been measured in S-type AGB stars and in meteoritic Stardust SiC grains 
that originated from AGB stars of roughly solar metallicity. The data of ^^Zr/^^Zr from 



meteoritic grains are more precise (with uncertainties as low as 5%— 10%, iLugaro et al. 



2003M) than the data 



limits ( iLambert et al. 



Tom S -type stars, where often it was only possible to derive upper 



19951 ). The differences of 22% found here when changing the 



^■^C(a;, n)^^0 rate do not drastically change the interpretation of the data, which mostly 
indicate strong deficits in ^^Zr, with respect to solar, and are overall matched by AGB 
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models (jLugaro et al. 



2003bD . In the case of K r/ Kr, Stardust SiC data range from ~0.5 



to 2 of the solar ratio of 1.52 (ILewis et al. 



19941 ). Our AGB models of metallicity around 



solar (Z = 0.01, 0.02), which are believed to well represent the parent stars of this S tardust, 



Verchovsky et al. 



(|200J) showed 



only produce ^^Kr/*^Kr ratios lower than solar. However, 
that the ^^Kr atoms are probably implanted in SiC by the high-velocity winds experienced 
during the evolution that follows the AGB, the post-AGB and the planetary nebula nucleus 
phases. The effect of the ^^C(a, r;,)^^0 reaction rate needs to be evaluated with regards to 
the composition of Kr inside the He intershell, which is exposed to the surface during these 
final phases of the evolution instead of the surface abundances as done in the present work. 
In this context the possibility of a proton ingestion in the very late TP sometimes occurring 
during the post-AGB evol ution also has t o be taken into account. First simulations of this 



proton ingestion event by 



Herwig et al. 



(120111 ) have shown that the ^^C{a,ny^O reaction 
rate plays a main role in determining the final abundances (see their Fig. 12). These models 
need to be further investigated in relation to the present ^^C{a, nY^O reaction rate and the 



86 



Kr/^^Kr ratios in SiC grains. 



5. Summary and conclusion 

We determined the stellar rate of the ^^C(a, r;,)^^0 reaction and incorporated the 
new reaction rate in calculations of the s-process nucleosynthesis in AGB stars. The 
Sa and ANC for the 6.356 MeV 1/2+ subthreshold state in ^^O were obtained from the 
measurement of the ^•^C(^^B, ^Li)^^0 angular distribution. This provided an independent 
examination to shed some light on the existing discrepancies in the and ANC values 
derived from different authors. Based on the measured ANC, we extracted the a-width of 
the 1/2"*" state in ^^O, which is currently the most uncertain parameter for determining the 
^'^C(a, ?7.)^^0 reaction rate. By using the present a-width and considering the properties 
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of ^^O states up to 8.342 MeV as well as their interferences, we derived the astrophysical 
S-factor and the stellar rate of the ^^C{a,ny^O reaction. At a temperature of 100 MK 
the new rate is roughly two times larger than the CF88 value and two times smaller 
than that recommended by NACRE (see Fig. [5]). Verification of the prese nt result using 



other independent techniques is desir able, e.g., the Trojan horse approach ( iSpitaleri et al 



1999 



Mukhamedzhanov et al. 



20081 ). and the iso spin symmetry approac h based on a 



20071 ). In addition, an 



measurement of the 1/2"*" 6.560 MeV state in ^^F (ITimofeyuk et al. 
extension of the experimental data of ^^C{a, riY^O toward lower energies is highly desirable, 
which can probably on ly be performed in an underground laboratory, e.g., LUNA (see 



Costantini et al. 



20091 ). 



We incorporated different ^^C(a,n)^^0 reaction rates in calculations of the s-process 
nucleosynthesis in AGB stars and found that: (1) If ^'^C burns completely in radiative 
conditions during the interpulse phase (as for stars of initial mass greater than ~2 Mq), 
there is no change in the final results. (2) If some ^'^C burns instead inside the convective 
TPs (for stars of initial mass lower than ~2 Mq), we find changes of up to 25% in the 
s-process results, particularly for Pb. There are model uncertainties related to the result 
of point (2): a) when ^^C burning in the TPs is due to incomplete burning of the ^'^C 
during the interpulse period, the exact stellar mass and metallicity range where incomplete 
burning of the ^^C during the interpulse period occurs, as well as TP numbers, and the 
exact amount of ^'^C ingested, all depend sensitively on the temperature and density in the 
^•^C pocket, on the interpulse period, and on the details of the inclusion of the ^'^C pocket, 
b) When ^'^C burning in the TPs is due to ingestion of protons directly inside the TP (as 
in l ow-mass and low-me tallicity stars, as well as in post-AGB stars experiencing a late TP, 



see 



Herwig et al. 



20111 ). the amount of ^'^C and neutrons produced strongly depends on 



the physical and numerical treatment of the mixing scheme adopted, which is at present 
uncertain. 
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Due to these model uncertainties together with the fact that the stellar observations 
have relatively large error bars, it is presently not possible to conclude if the new rate 
of the ^^C{a,ny^O reaction provides the best match to the available observational 
constraints. This may be however possible in the future, when development of recent 3D 



201 If ) will allow 



hydrodynamical models of the proton ingestion episodes ( iStancliffe et al. 
a better understanding of neutron production and the s-process inside TPs to be compared 
to the composition of stellar observations and Stardust grains. 
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Fig. 1. — Focal-plane position spectra of the ^Li events at ^lab = 10° from the a-transfer 
reactions. The solid and dashed lines denote the results from the enriched ^^C target and 
natural ^^C target, respectively. 
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Fig. 2. — Angular distributions of the ^^B+^^C elastic scattering at incident energy of 50 
MeV and the '''Li+^^0 elastic scattering at incident energy of 26 MeV. The solid curves 
represent the calculations with the fitted OMP parameters. 
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Fig. 3. — Angular distributions of the ^^C{^^B, ^Liy^O reaction leading to the excited states 
at = 3.055 MeV, 3.843 MeV, 4.554 MeV and 6.356 MeV. The curves denote the DWBA 
calculations with the fitted OMP parameters. DWBAl and DWBA2 repr esent the results 
using the spectroscopic amplitudes of ^^B from 
respectively. 
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(120081 ). respectively. The solid and dotted 



(red) curves represent the present S-factor and its uncertainties. The dashed (blue) curve 
represents the S-factor by excluding the contribution of the 1/2"*" subthreshold resonance. 
(A color version of this figure is available in the online journal.) 
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of 0.1 GK. (A color version of this figure is available in the online journal.) 



- 34 - 




time 



Fig. 6. — Schematic diagram of the evolution of the structure of an AGB star with the 
superimposed locations where neutrons are released by the ^"^C neutron source in the different 
regimes described in the text. Also indicated are the masses of the models and the typical 
TPs where each regime occurs. (A color version of this figure is available in the online 
journal.) 
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Table 1: OMP parameters used in the present DWBA calculation. E^^ denotes the incident 
energy in MeV for the relevant channels, V and W are the depths (in MeV) of the real and 
imaginary parts of Woods-Saxon potential, and r and a are the radius and diffuseness (in 
fm) of Woods-Saxon potential. 

Channel E^^^ V W 
i^B+i^C 50.0 182.64 0.788 0.740 8.193 1.250 0.740 
^Li+i^O 26.0 114.20 0.737 0.719 34.602 0.997 0.764 
^Li+i^C 31.8 159.00 0.630 0.810 8.160 1.330 0.780 



Table 2: The Sa and ANC (C^) as derived from different measurements for the 1/2^*" sub- 
threshold state of ^^O. 



Reference 


Transfer system 


Incident energy 




^2 






(MeV) 




(fm-i) 


KubonoetaL r2003) 


i3C(6Li,rf) 


60 


0.011 




Keelev et al. f2003) 


i3C(6Li,rf)^ 


60 


0.36-0.40"^ 




Johnson et al. r2006) 


6Li(i3c,rf) 


8.0 and 8.5 




0.89 ± 0.23 


Pellemti et al. r2008) 


i3C(^Li,t) 


28 and 34 


0.29 ± 0.11 


4.5 ± 2.2 


Present work 


i3C("B,7Li) 


50 


0.37 ± 0.12 


4.0 ± 1.1 



"Reanalysis of the experimental data of iKubono et al. I (|2003l l. 

^See Table 3 of iKeelev et al. I (|2003l ). these values also depend on different normalization procedures used in 
their work. 
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Table 3: Resonant parameters for ^''O used in the present calculations. and Ej^ represent 
level energy and resona nt energy, respective ly. E^;, Er, J'", Tn{Eji) and Ta{Eji) are taken 
from the compilation by ixilley et al. ( 1993 ). 



Level 


E, (keV) 


Er (keV) 




Tn{ER) (keV) 


To,{Er) (keV) 


7a (keV) 


1 


6356 


-3 


1/2+ 


124 




12. 7^^ 


2 


7165 


806 


5/2- 


1.38 


0.0033 




3 


7201 


842 


3/2+ 


340*^ 


0.08^ 




4 


7378 


1019 


5/2+ 


0.64 


0.01 




5 


7381 


1022 


5/2- 


0.96 


0.003 




6 


7558 


1199 


3/2- 


500 


0.08 




7 


7688 


1329 


7/2- 


13 


0.01 




8 


7956 


1597 


1/2+ 


84 


6.7 




9 


7991 


1632 


1/2- 


250 


16 




10 


8058 


1699 


3/2+ 


71 


15 




11 


8200 


1841 


3/2- 


48 


4.0 




12 


8342 


1983 


1/2+ 


8.1 


2.2 





"Reduced a-width from this work. 

''The a- and neutron- widths for the 3/2+ resonance at Eji — 842 keV were adjusted to provide the best fitting 



of the experimental data, which resuhed in broader 
280 keV, TaiER) = 0.07 keV] of 



Tillev et al. 



partial widths than the recommended values {TniEji) = 



( 19931) . 



- 38 - 



Table 4: Numerical values of the present ^^C{a,nY^O rates (cm^/mol/s) with the adopted 
value, upper and lower limits for the 0.04—3.0 GK temperature range and the coefficients Oj 
in Eq. [8] for these three rates. The overall fitting errors are all less than 7% at temperatures 
from 0.04 to 10.0 GK. 



n 




Reaction rate 










Coefficients 




Adopt 


Upper 


Lower 




ai 


Adopt 


Upper 


Lower 


0.04 


2.63x10-24 


3.41xl0-2'* 


1.92x10- 


24 


ai 


79.3008 


52.8016 


87.5453 


0.05 


1.63x10-21 


2.07x10-21 


1.23x10- 


21 


a2 


-0.304890 


-0.204024 


-0.301248 


0.06 


2.17x10-19 


2.70x10-19 


1.67x10- 


19 


as 


7.43132 


-23.2591 


9.17396 


0.07 


1.07x10-1^ 


1.31x10-1^ 


8.44x10- 


-18 




-84.8689 


-42.4710 


-95.7947 


0.08 


2.67x10-16 


3.21x10-16 


2.14x10- 


16 


ag 


3.65083 


35.1371 


4.66751 


0.09 


4.04x10-1^ 


4.78x10-15 


3.28x10- 


-15 


as 


-0.148015 


-17.3173 


-0.221941 


0.10 


4.19x10-1^ 


4.90x10-1^ 


3.46x10" 


-14 


ay 


37.6008 


6.45708 


40.9578 


0.15 


1.67x10-1° 


1.86x10-1° 


1.46x10- 


10 


as 


62.5775 


64.3536 


63.1390 


0.20 


3.33x10-^ 


3.59x10-^ 


3.03x10- 


-8 


ag 


-0.0277331 


-0.302435 


0.0195985 


0.30 


3.09x10-5 


3.20x10-5 


2.96x10- 


-5 


aiQ 


-32.3917 


3.36966 


-34.5026 


0.50 


8.06x10-2 


8.12x10-2 


7.98x10- 


-2 


an 


-48.9340 


-64.0633 


-47.2196 


1.0 


2.44x102 


2.45x102 


2.44x10 


2 


ai2 


44.1843 


1.62313 


44.0189 


2.0 


5.99x10^ 


6.00x10^ 


5.97x10 


4 


an 


-20.8743 


0.00566612 


-20.9558 


3.0 


6.12x10^ 


6.13x105 


6.10x10 


5 


0-14 


2.02494 


31.1730 


0.905521 
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Table 5: Summar y of the stellar models. 

Mass (Mq) Metallicity (Z) Number of ^^C pockets 
1.0 0.0001 2 

1.5 0.0001 15 

1.8 0.01 6 

3.0 0.02 16 



